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0!7  Till.  MUATI0H3  CV  HUt-lRICAL  PREDICTION 


IN  RELATION  TO  THE  KON-DIVSROSNCE  CONCEPT 


by 

V.V.  Bykov 


The  numerical  method  of  short  range  weather  prediction 
makes  use  of  a  system  of  equations  consisting  of  two  equations, 
of  motion,  the  equation  of  continuity,  the  equation  of  heatflow 
(in  adiabatic  conditions)  and  the  hydrostatic  equation.  Taking 


the  horizontal  co-ordinates  x  and  y,  the  reduced  pressure 

whore  p  is  the  pressure  and' P  the  pressure  at  sea  level) 
t-tho  time,  as  independent  variables,  the  above-mentioned 
equations  nay  be  written  as  follows: 


and 


Q 

^  Rer\ 


(1) 

& +«£+•$— 

(2> 

iu  ,  dv  ,  1  dt 

■$rr?j+p?r=0' 

(3> 

*T  A-..dT  -L.„dT  *  O. 

«> 

T  C  9H 
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and  v  - 
H  - 
T  - 
1  - 
R  - 
W  *£.  . 

.  at 


wind  components  along  the  x  and  y  axis 

contour  height 

temperature 

Coriolis  parameter 

gas  constant 

Individual  change  of  pressure,  playing  in  our 
system  of  independent  variables  tho  r ole  of 
vertical  motion,  being  related  to  it  by 


where 


w  -  the  vertical  component  of  wind 
g  -  the  acceleration  of  gravity 


C'—Rt, 


where  T«-=  the  dry  adiabatic  lapsrate 

T  =  tho  mean  lapsrate 

f  s  the  mean  temperature 

— N  The  system  of  equations  (l)-(5)  is  solved  with 

respect  to  the  initial  and  boundary  conditions.  Tho  "t  and  T 
can  oa3ily  be  eliminatod  from  eq.(3)-(5)  and  the  resulting 
equation  systei:  consists  of  throe  equations  with  throo  unknown 


functions  u,v,  and  H.  The  latter  equation  system  is  solved  for 
geostrophic  conditions,  i.e.  it  is  assumed  that  horizontal 
motion  ot  atmospheric  disturbances  and  the  temperature  field 
takes  place  geostrophically ,  tho  wind  components  u  >=>  and 

l  dH  f  9 

VmtY'Sx  *>eing  determined  uniquoly  by  the  contour  field.  For 

that  reason  the  latter  is  taken  as  the  initial  state.  The 


boundary  conditions  are  usually  taken  to  be  that  at 


and 


at  C— 1  t 


•PoH 
t*  IT  *' 


C -t>* 
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The  validity  of  the  application  of  the  geostrophic 
approximation  is  demonstrated  in  many  papers  (3,4,5).  In 
particular,  I. A.  Kibel  (4)  has  given  a  generalized  method  for 
the  adaptation  of  atmospheric  motion  to  the  geostrophic  motion. 


However,  the  general  acceptance  of  the  geostrophic 
approximation  and.  its  application  in  the  equations  of  numerical 
prediction  has  not  meant  the  end  of  studies  on  the  doparturosof 
the  actual  wind  from  the  geostrophic.  Many  papers  (6,11,12)  deal 
with  the  estimation  of  errors  ;rosulting  from  tho  application  of 
the  geostrophic  concept  in  the  equations  of  atmospheric  dynamics. 
In  some  equations  of  numerical  prediction,  the  non-goostrophlc 
departures  of  wind  have  been  accounted  for  by  their  approximate 
values  in  terms  of  geostrdphic  wind. 


An  alternative  approach  for  talcing  into  account  tho 
non-goostrophic  effect  in  the  :s,tudies  on  atmospheric  processes 
is  to  accept  the  hypothesis  of  non-divergence  of  flow,  based  on 
the  smallness  of  the  horizontal'  velocity  divergence.  In  this, 
scheme  it  is  assumed  that  the  wind  field  can  bo  expressed 
approximately  et  r  11  JL^V':l3  of .  the  atmosphere  by  a  stream 
function^;,  so  thatB^_*fc  •  d*  .This  approach  has  beon  suggested 
in  (1)  arid  (2).  ST*- 


In  recent  years  3olin  (?)  and  Charney  (0),  on  tho  basis 
of  oxporlonce  in  numerical  prediction  with  electronic  computers, 
have  suggested  that  tho  concept  of  non-divorgonce  be  used  in  the 
equations  of  numerical  prediction  instead  of  geostrophio  approxi¬ 
mation  as  a  means  of  inoreasing  computational  stability  in  fore¬ 
cast  ing . 

Accepting  this  concept  and  introducing  the  stream 
function  into  the  equations  of  numerical  prediction,  thoso  latter 
can  written  for  tho  contour  field  II  at  the  level  of  non-divorg- 
enco’  us  follows: 


where  p  “J^T* 

This  system  is  solved  for  a  limited  flat  region  using 
for  tha  initial  state  the  H(x,y)  values  from  the  contour  map 
that  corresponds  to  the  ievel  of  non-divergenoe . 


Tho  method  of  solution  of  those  equations  is  as  followo: 
Tho  non-linear  equation  (7)  with  known  rlght-hond  side  is  first 
solved  for  •  This  yields  the  initial  stream  function  correspond¬ 
ing  to  the  initial  H  field.  The  solution  of  oq.(7)  is  then 
introduced  Into  eq.(6)  and  this  la  aolved  for‘$  ,  yielding  tho 
future  distribution  of  the  stream  function.  This  is  put  back 

*  According  to  LanderB  (10)  this  level  is  beat  represented  by 
the  700  mb  surface. 
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again  into  oq.{7)  and  the  latter  Is  solved  for  the  now  unknown 
function  H,  i.e.  the  future  contour  field  is  calculated  from 
the  future  stream  function. 


The  basic  difficulty  in  the  practical  application  of 
this  model  is  the  initial  solution  of  eq.(7),-  the  equation  of 
balance  being  a  non-linear  differential  eq'uatlon  of  Monge- 
Ampere  type. 


Let  us  examine  the  individual  terms  in  this  equation. 
If  the  non-linear  terms  and  the  torm  for  the  variation  of  the 
Coriolis  force  with  latitude  are  neglected,  then  the  term 


7  dtf,  relates  the  oontour  field  with  the  stream  function  in 
a  quasi-geostrophic  approximation.  In  other  words  the  eq.(7) 
transforms,  on  the  given  assumption,  into  the  equation  of 

vortiolty:  q _ £*i . 

dx  dy* 

cK» 

The  tern  jPjy'  in  eq.(7)  indicates  that  the  application 
of  the  goostrophic  approximation  for  vortiolty  results  in  an 
underestimate  of  vortiolty  in  westerly  flow  (^<0)  and  ftn 

overestimate  in  the  easterly  flow  (^>0j.( 

The  non-linear  terms  “  (3737 )*  can  have  a  pronounoed 

of foot  only  in  oase  of  a  large  ourvature .combined-  with  largo 
horizontal  gradients,  as  observed  in  well  developed  oyclonesand 
anticyclones  and,  to  a  lesser  extent,  in  troughs  and  ridges. 

In  particular ,  the  neglect  of  the  non-linear  terms  (whioh  occurs 
with  the  quasi-geostrophio  approximation)  leads  to  an  ovor- 
ootimato  of  vortiolty  in  oyolones  and  an  underestimate  in' anti¬ 
cyclones.  Some  idea  of  the  aotual  magnitude  of  these  terms  of 
the  Monge-Ampere  equation  oan  be  obtained  from  Fig.l. 


Lot  us  now  disouss  some  peculiarities  of  the  numerical 
solution  of  eq.(7).  This  equation  Is  of  the  elliptical  typo  with 
rospeot  to  the  unknown  funotion' 4"( x,y)  ,  if  the  condition 


**- ?£>— r- 


is  satisfied  everywhere  in'  and  on  the  boundaries  of  the  region 
undor  consideration. 


This  condition  is  fulfilled  practically  everywhere 
on  oho  particular  contour  map  that  corresponds  to  the  level  of 
non-divergence .  For  the  numerical  solution  of  oq.(7),  Bolin  (7) 
suggests  tho  relaxation  method  for  one  variant  of  consecutive 
approximation.  The  speed  of  convergence  of  the  iteration 
process  In  solving  a  differential  equation  by  relaxation  depends, 
as  ia  known,  on  how  successfully  we  have  seleoted  the  first 
approximate  solution.  The. .latter  for  eq.(7)  oan  be  taken,  for 
instance,  in  the  f orm  «[>  =  y,  i.e.  it  is  assumed  that  the  stream 
function  is  related  to  the  contour  field  as  in  the  case  of 
geootrophio  approximation. 


The  second  variant  for  approximate  solution  of  eq.(7) 
night  be  obtained,  for  inotanoe,  in  the  following  way.  Lot  us 
rewrite  eq.(7)  as  follows: 


w.  „  lot  uo  consider  it  as  a  Poisson  equation  for  tho  unknown 
function  <7 *=  <}  — -y-  with  known  right  hand  side,  which  expresses 
the  vorticity  difference  assuming  a  non-divergont  flow  instead 
of  a  geostrophlo  flow. 
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Eq.(9)  is  solved  numerically  in  the  follo\»ing  way. 

J'/e  rewrite  it  In  finite  differences  for  n  gridpolnts,  symmet¬ 
rical  with  respect  to  a  central  point  (i,j),  for  which  the 
solution  is  sought.  On  the  boundaries  of  the  region  under  con¬ 
sideration  we  assume  the  flow  to:  be  geostrophic  jl.  e .  \g  ^ _ it  s  o 

In  this  way  one  obtains  a  system  of  n  linear  algebraic  equations 
with  n  unknowns.  The  solution  of  the  above  system  for  the  unknown 
q  at  point  (l,j)  represents  a  sum  of  products  of  the  right  hand 
side  of  eq.(9)  £o,(2C,y)J,  taken  at  the  surrounding  n  points 
(including  also  the  oentral  point)  by  the  corresponding  weight 
fuctors.  The  right  hand  side  of  eq.(9)  is  computed  on  the 
assumption  of  the  geostrophic  wind,  i.e.  jo  — -y-.j 

The  values  of  the  individual  terms  and  the  whole  right 
hand  side  of  eq.{9),  computed  from  the  AT7OO  map  for  the  morning 
of  Sept.  7,  1953 ,  are  given, in  Fig.  1.* 

The  method  outlined  here  for  the  short  range  numerical 
prediction  of  the  contour  field  at  the  level  of  non-divergonco, 
based  on  the  solution  of  eq.(6)  and  (7)  assuming  non-divergent 
flow,  can  be  realized  in  prnotioe  by  means  of  electronic  computers 
as  follows: 


The  initialitfa'  at  each  grid  point  is  read  off  from 
an  analyzed  oontour  map' for  the  level  of  non-divergence.  Those 
are  inserted  into  eq,(7)  and.,  the  latter  is  solved  forjy/,  yielding 
the  initial  stream  function^;.  This  is  now  introduced  into  eq. 
(6)  which  forooasts  the  future  stream  funotion: 


« 


iK 

~Sy' 


The  actual  computation  of  the  future  stream  function 
is  porformed  by  time  steps.  Introducing  into  eq.(6),l*Lj 

and  the  whole  right  hand  side,  containing  derivatives  of  tyYwith 
rospoot  to  the  coordinates,  is  computed  first  for  the  timo  instant 
i.e.  the  quantity  j~l  is  calculated  for  the  time  instant  7# , 

assuming!  ^  —  6  •  and  —  at  the  boundaries  of  the  forecast  region. 


After  computation  of 


<**>«•» 1 -<*«<, +  (t?) 


the  equation 


(10) 


is  applied  at  all  gridpolnts,  where  (M1  is  the  "time  step"  selected, 
usually  according  to  the  dimensions  of  the  grid. 


After  the  numerical  v  alue  of;dij<,'at  the  Instant  (/er+*V>) 
is  computed  for  the  internal  grid  points,  i.e.  tho  right  hand 
Bide  of  eq.(10)  is  numerically  solved  for  ^  as  a  Poisson  equation, 
i.e.,  in  tho  form  ( ^  (*»  £)»!.  the  :iew  value  of  stream  function 
W/.t-a  itoelf  is  calculated  for  the  grid  points. 

The  A<jP;  and  .<(>;  obtained  for  time  instant  ty-h*.  are 
now  lnsorted  into  oq.(10)  and  the  analogous  procedure  is 
repeated  for  the  time  instant  2  dt  etc,  etc. 


After  the  (&)/.+«/  and :  (<*'!>), are  found  for  a  given 
time  interval  [kit,  where  k  is  tho  number  of  time  3tepe,  they 
are  introduced  into  eq.(7). 


Considering  this  equation  now  as  a  Poisson  equation 
for  the  unknown  function 


#  This  figuro  is  not  reproduced  in  this  translation. 
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with  known  right  hand  side,  it  is  solved  numerically  for  H, 
i.e.  the  future  contour  field  is  calculated. 

/  Next,  we  derive  the  system  of  prognostic  equations 

for  a  generally  baroolinic  atmosphere  assuming  a  non-divergent 
windfield.  Nor  this  wo  make  use  of  the  initial  equations  (l)-(S/. 

Differentiating  eq.(l)  with  respect  to  y  and  eq.(2) 
with  respect  to  x,  and  subtracting  the  first  from  the  second,  one 
obtains,  as  usual,  the  vorticity  equation.  Applying  now  eq.(S) 
for  continuity,  we  write  the  resulting  equation  in  tho  following 
form: 

£9  +  uia  +  *i9.  +  8t,_4.±L-0  (12) 

dt'  ox  ~  T r  P  <K  u  * 


whore  C  — -jy is  the  vertical  component  of  vorticity,  and 

tho  X-axi 8  is  directed  toward  the  east,  tho  y  axis  toward  tho  north 
In  addition,  wo  have  disregarded  (1  in  eq.(12)  as  negligibly 
small  in  comparison  with  1.  . 

Differentiating  eq.(l)  with  respect  to  x  and  eq.(2) 
in  rospect  to  y  and  summing  up,  one  obtains  the  equation  of 
divergence: 

^  +  +  +  +  —  .  U3) 


, where  P-mm'7rz~r'j%‘  i°  the  horizontal  velocity  divorgonce. 

No  restrictive  assumptions  havo  so  far  been  made  in  respect  of 
wind  speed  in  oq.(12)  and  (13). 

Let  us  further  assume,  that  the  condition 

1141 


is  fulfilled  approximately  at  all  levels  of  the  atmosphere.  The 
horizontal  velocity  divergence  is  thus  assumed  to  bo  approximately 
zero  everywhere,  exoept  for  the  vortiolty .  oquation  ( term  _ 

This  assumption  is  based  on  the  estimates  of  the  oraor  of  magni¬ 
tude  of  terms  in  the  hydrodynamic  equations  for  tho  atmosphere. 

The  non-divergence  (quasi-solenoidal)  assumption  for 
ntmoopheric  motions  (eq.14)  allows  one  to  introduce  a  stream 
function,  so  that 

“  “  ~77>  VW3-33T-  (15) 

■  n  dv  da  . , 

than,  i  “***57  — and  assuming  the  wind  field  is  non-divergent, 
oq . ( 12 )  and  eq.(l3)  may  be  written  as  follows: 


o _ dv  da.  . . 

°*"3 7-37 


6 


In  addition,  we  transform  eq.(4)  for  heat  flow.  Using 
eq .  (5)  and  (15)  we  write  it  in  the  following  form: 

*[*($+(*  S)!l+*'-o-  1101 

Differentiating  eq.(18)  with  respeot  to  C 

(15> 

and  eliminating  by  combining  equations  (16)  and  (19)  we  have: 

’•  *  •  ,  ,  % 

Eq.(20)  oontains  two  unknown  functions  and 

The  seoond  equation  relating  these  unknowns  is  obtained 
by  differentiating  eq.(17)  with  respect  to  time  t.  The  resulting 
equation: 

\  *(4W[3$(»+8£(*)-  U1) 


has  variable  ooeffioientB,  however,  it  is  like  eq.  (20)  linear  with 
respoot  to  unknown  and 

Thus,  the  system  of  eq.(20)  and  (21)  can  be  used  for  the 

O determination  of  unknown  fvinotions  -and  j3L  In  addition  to  eq.(2C) 

7T  dt' 

and  eq. ( 21) ,  eq.(l7)  is  used,  in  analogy  with  the  previous  modol 
for  the  level  of  non-dlvergenoe,  to  determine  the  initial  ^ 
funotlon  from  the  given  H  field.  Two  boundary  conditions  for’C. 
are  given  in  the  above  equation  system. 

The  solution  of  the  equation  system. eq.( 17) ,  (20)  and  (21), 
l.o.  the  determination  of  the  future  distribution  of  wind  field 
(stream  funotlon)  and  contour  field  oan  bo  realized  in  practice  by 
means  of  eleotronlo  computers. 


r 
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